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ABSTRACT 



We propose a novel method for the construction of orthogonal arrays. The algorithm 
makes use of the Kronecker Product operator in association with unit column vectors 
to generate new orthogonal arrays from existing orthogonal arrays. The effectiveness 
of the proposed algorithm lies in the fact that it works well with any linear seed 
orthogonal array without imposing any constraints on the strength or the number of 
levels. The resulting orthogonal array has the same strength as the seed orthogonal 
array. We also discuss the proof of correctness of the algorithm. In the Results section 
we provide a list of new orthogonal arrays generated using this algorithm, that are 
currently not present in the libraries of orthogonal arrays. 
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Chapter 1 
Introduction 



Construction of orthogonal arrays (OAs) is an important problem in combinatorial 
design which holds great significance for design of experiments in statistical analy- 
sis. In the past, several construction methods for generating orthogonal arrays have 
been proposed and analyzed. Hedayat et al. (1999) provide a comprehensive study 
of orthogonal arrays. The significance of the Kronecker Product and the Kronecker 
Sum operations in the context of generating orthogonal arrays is well established, as 
seen in the works of Shrikhande (1964), Wang and Wu (1991), Zhang, Weiguo, Mao 
and Zheng (2006) and Sinha, Vellaisamy and Sinha (2008). Sinha et al. (2009) used 
the Kronecker Sum operation on ternary orthogonal arrays and Balanced Incomplete 
Block Designs to construct new symmetrical ternary orthogonal arrays. In the cur- 
rent work we propose a novel construction approach for orthogonal arrays using unit 
column vectors and the Kronecker Product operations on existing orthogonal arrays. 
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1.1 Preliminaries and Definitions 



1.1.1 Balanced Arrays 

A balanced array, denoted by BA{N, m, s, t) {lJ>xi,....,xt}j is defined as an A?" x m matrix 
B with elements belonging to the set S = {0, 1, s — 1} of s symbols, m factors, N 
runs and strength t such that every N x t sub-matrix of B contains the ordered row 
vector {xi, ....,Xt), /J-x^.-.^xt times, where Hxi,....,xt is invariant under any permutation 

of Xi, .....Xf 

1.1.2 Orthogonal Arrays 

An N X k array A with entries from S is said to be an orthogonal array with s levels, 
strength t {0 < t < k) and index A if every N xt sub-array of A contains each i-tuple 
based on S exactly A times as a row. 

H lJ'xi,....,xt — A*(constant) V i — tuples {xi, . . . ,Xt) G 5" in a balanced array, then 
the balanced array becomes an orthogonal array with index /i. 

1.1.3 Kronecker Sum 

Sinha et al. (1979) defined the Kronecker sum of matrices A (order m x n) and B 
(order p x q) asA^B — A^J+J^B , where (g) denotes the usual Kronecker 
product and J is a matrix with all its elements unity but is of dimension p x ^ in the 
first term and m x n in the second term. The Kronecker sum of binary orthogonal 
arrays was defined by Sinha et al. (2008). 

1.1.4 Galois Field 

A Galois field is a field that contains finite number of elements and is denoted by 
GF{s), where s is its order. The order of a Galois field is the number of elements in 
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the field and is of the form p", where p is a prime number known as characteristic of 

the field and n is a positive integer. We shall denote the elements of the GF{s) by 
{0, 1, 2, s — 1} and the set of all n-tuples with entries from GF{s) by GF(s)". 

1.1.5 Simple Orthogonal Arrays 

An orthogonal array is said to be simple if all its runs are distinct. 

1.1.6 Linear Orthogonal Arrays 

Let s = p", where p is a prime and n is a positive integer. Then, the orthogonal 
array OA{N, k, s,t) with levels from GF{s) is linear if it satisfies the following two 
conditions: 

• it is simple 

• when its rows are considered as /c-tuples from the GF{s), its N runs form a 
vector space over GF{s). 
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Chapter 2 



Related Work 



The problem of generating orthogonal arrays has been of interest to researchers since 
more than half a century. A number of algorithms have been proposed and construc- 
tions given. The orthogonal arrays that have been successfully constructed using such 
constructions are stored in libraries. In this chapter we shall discuss some important 
results in this area of research so as to build a basic framework before discussing the 
work presented in this thesis. We shall discuss the theorems and constructions briefly 
without digressing by going into the details of the proofs. For detailed proofs the 
reader is directed to the references. 



2.1 The Rao-Hamming Construction 

As the name indicates, this construction was given by Rao (1947, 1949) and Hamming 
(1950), both of whom had found this algorithm independently. Rao had originally 
introduced a special case of the idea of orthogonal arrays in a rather implicit sense in 
his concepts of hypercube of strength t. His construction of hypercubes of strength 2 
find relevance and correspondence with orthogonal arrays OA{s"', (s** — l)/(s— 1), s, 2), 
s being a prime power. This brings us to the following theorem. 
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Theorem 2.1: If s is a prime power then an OA{s'^, {s'" — — 1), s, 2) exists 
whenever n > 2. 

Construction: Let us consider an s" x n array whose rows are all possible n- 
tuples over GF{s). Let, Ci, C2, , Cn be the columns of this array. The columns of 

the orthogonal array then consist of all the columns of the form 

ZlCi + Z2C2 + + ZnCn — [Ci, C2, , Cn]z, 

where z = (^1,^2, > ^n)^ is an n-tuple from GF{s), not all the Zi are 0, and the 

first 2;^ is 1. Then, there are (s"" — l)/(s — 1) such columns. 

2.2 Bush's Construction 

Bush's (1952) research on orthogonal arrays of index 1 is a well known and important 
result. The theorem proposed by Bush goes as follows. 

Theorem 2.2: If s > 2 is a prime power then an OA{s*, s + l,s,t) of index unity 
exists whenever s > t — 1 > 0. 

The result stated in the above theorem can be improved for some values of s and 
t as can be seen below. 

Theorem 2.3: If s = 2'",m > 1, and t = 3 then there exists an OA{s^, s + 2, s,t). 

The orthogonal arrays constructed in Theorems 2.2 and 2.3 are simple and hnear. 

Theorem 2.4: If s is a prime power and a hnear array k, s, t) exists, then 

there also exists a linear array OA{s'^~*, k,s,k — t). 

2.3 Bose and Bush's Recursive Construction 

This construction was proposed by Bose and Bush (1952). It allows for the construc- 
tion of orthogonal arrays of strength two with a large number of factors and possibly 
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the the maximal number, provided that the number of symbols s and the index A are 
powers of the same prime. The theorem is stated as follows: 

Theorem 2.5: Let s — p" and A = p", where p is a prime and u and v are integers 
with u>Q,v>l. Let d— \u/v\. Then there exists an 



2.4 Hadamard Matrices and Orthogonal Arrays 

A Hadamard matrix is a square matrix which takes only two symbols +1 and -1 as its 
entries such that for every two different rows there are matching entries in exactly half 
of the cases and non-matching entries in the remaining half. It can be easily noticed 
that Hadamard matrices are difference schemes with two symbols. Hadamard matri- 
ces and orthogonal arrays have close resemblances in their combinatorial properties. 
Hadamard matrices can be generated using recurrence relations. One such method is 
called Sylvester's method which goes as follows: 

Let Hi denote a Hadamard matrix of order i. Then, 



d+l 



l)/{s'-s''-') + l,s,2) 



Hi 



1 



1 



1 



H2^ 



1 -1 




H2k 
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An important result that illustrates the close connection between Hadamard ma- 
trices and orthogonal arrays is mentioned below. 

Theorem 2.6: Orthogonal arrays 0A(4A, 4A - 1, 2, 2) and 0^(8A, 4A, 2, 3) exist 
if and only if there exists a Hadamard matrix of order 4A. 
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Chapter 3 

A Novel Construction Algorithm 



In this chapter, we propose a novel method for the construction of orthogonal arrays 
using existing orthogonal arrays and unit column vectors with the help of the Kro- 
necker Product operator. The method proposed herein serves towards forming orthog- 
onal arrays of larger dimensions from orthogonal arrays of smaller dimensions. We 
primarily deal with linear orthogonal arrays, however as we will show later, non-linear 
orthogonal arrays can also be constructed using this approach. We then investigate 
the correctness of the construction and other possible extensions to the algorithm. 

3.1 The Construction 

Let A denote a hnear seed OA{N,k,s,t) with levels from GF{s). The construc- 
tion provides an approach to generate another OA{N'^,k'^ + 2/c,s,t). Let C = 
{ci, C2, . . . , Cfc} denote the set of all the factors of A. Now let us define another 
zero column vector Ck+i- Let U denote a unit column vector of dimensions N x 1. 
See Algorithm 1 for the detailed construction. 

Algorithm 1 returns a set of column vectors C. Then, C forms the set of all the 
factors of an OA(N'^, k'^ + 2k, s, t). Let us denote this OA(N'^, k'^ + 2k, s, t) as B. 
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Algorithm 1: The Construction Procedure on (A,C,U,s) 



begin 

i ^ 1 

C = CU{ck+i} 

while i < A; + 1 do 
while j < k + 1 do 

ii i ^ k + 1 or j ^ k + 1 then 

c' = Cii^U + U <^Cj 
d — c'(mod s) 

a ^c'u {d} 

L i i + 1 
_ return C 



3.2 Correctness of the Algorithm 

Now let us discuss the proof of correctness of the construction proposed above. We 
shall state the result formally in the form of the following theorem. 



Theorem 3.1: The existence of a linear orthogonal array OA{N,k,s,t) with 
levels from GF{s) impUes the existence of an orthogonal array OA{N'^, k"^ + 2k, s, t). 



Proof. We proceed with the proof by first showing that the OA{N^, k^ + 2k, s, t) 
generated from the seed hnear orthogonal array OA{N, k, s, t) is itself linear. Recall 
the denotations used in the previous section since they will be used again in this 
proof. 

We can write B = [Bi\B2 \ . . . \Bk\Bk+i\Bk+2]-, where Bi is defined as follows: 



B, 



— {ci®U + U ® Cj)mod s, Vj = 1 to k} 
{Cfe+i|c'fc+i = ([/ (8) Cj)mod s, Mj = 1 to A;} 



{4 



k+2\'-'k+2 



{cj ® U)mod s, Vj = 1 to A;} 



for i = 1 to k 
if i = A; + 1 
if i = A; + 2 



Let Rj — [lyil^l ■ ■ ■ \K+2] be the row of B, such that 6^, 6^, ... , 6^_^2 
rows oi Bi, B2, ■ ■ ■ , Bk+2 respectively. Then due to the definition of the construction, 
it can be seen that the following results hold: 

• for i = 1 to fc, 6i = r + /3jt/ixfe 

• h+i = r 

where r and [/3i|/32| ■ ■ ■ |/?fe] are rows of A and C/ixik is a unit row vector. 

Thus we can find rows ri, r2, [/3i|/32| • • • and [/3i|/32| • • • \l^'k\ of A which give the 
the i*'* and j^^ rows of B as follows: 

Ri = . . . = [n + ^iC/ixfe|ri + ^2C/ixfe| ■ ■ ■ |ri + ^ikC/ixfe|ri 1^1/^21 ■ ■ ■ \Pk\ 

Rj = [b>\bi\ . . . H+2] = h + P[U,^k\r2 + P2Uixk\ ...\r2 + /^'^U^Mmi ■ ■ ■ W'kl 

A linear combination of the above two rows using scalars hi and /i2 from the 
GF{s) will result in the following: 

hiRi+h2Rj = [hiri+h2r2+(hipi+h2p[)Uiy,k\hiri+h2r2+(hip2+h2p'2)Uiy,k\ ■ ■ ■ \hiri+ 
h2r2 + {hxh + h2/3'k)Uiy,k\hin + /i2r2|/ii/3i + h2(3[\hil32 + /i2/32l • • • \hif3k + Mfc] 

Now, since A is known to be linear, and ri and r2 are rows of A, hiri + /i2'"2 is 
also a row of A. Similarly, [hi(3i + h2/3[\hi/32 + /i2/32l • • • l^i/^fc + ^2/3fc] is a row of A. 
Hence, by definition of the construction and its properties discussed earlier in this 
proof, hiRi + h2Rj is a row of B. It can also be easily noted that one row of S is a 
zero vector of 1 x (A:^ + 2k) size. Hence, it can be concluded that the rows of B form 
a linear subspace of GF{s)''. 

So, let N"^ — and let G be the generator matrix for B of dimensions n x (k'^+2k) 

such that the rows of B consist of all /c-tuples riG, where r] G GF{s)"'. Now, suppose 
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there exist t columns of B that are hnearly dependent over GF{s). Then, due to the 
nature of the construction, there exist x {1 < x < t) columns in A that are linearly 
dependent, which is a contradiction. Hence, every selection of t columns of B are 
linearly independent over GF(s). So, let us choose t columns of B and let Gi be 
the corresponding submatrix of generator matrix G. Then, the columns of Gi will 
be linearly independent. Also, the number of times a i-tuple r is present as a row 
in these t columns of B is determined by and is equal to the number of rj such that 
vGi = T. 

Now, since Gi has rank t, the number of such r] is for all r. Hence, B is an 

orthogonal array of strength t. This concludes the proof. 

3.3 Binary Orthogonal Arrays 

In this section we propose the following lemma, wherein we show that the linearity 
condition in Theorem 1 can be dropped for binary seed orthogonal arrays (i.e. OAs 
of strength 2) . 

Lemma 3.1: The existence of an OA{N,k, s,2) implies the existence of an 
OA{N'^,k^ + 2k,s,2). 

Proof. Let A be OA{N, k, s, 2) as described before and let c[ and C2 be any two 
distinct factors of B. Also, let us assume c[ and were generated from the factors 
Cj^, Cjj, Cj2, of A using the following equations: 

c'l — {ci^ <S> U + U <^ Cj^)mod s 

C2 = (cjj ®U -\-U ® Cj^)mod s 

Let c'j and C2 be denoted by [Xi : X2 : . . . : X^]' and \Y-i -.¥2 : . . . : Yjq]' , where the 
dimensions of X^ and 1^, Vi = 1 to n, is x 1. 
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In order to show the vahdity of the algorithm we spht the analysis into four cases 
as follows. 

Case 1 : Cj^ ^ cj^ such that Cj^, cj^ ^ Ck+i 

We know that Xi — (cj^ + piU)mod s and Yi — {cj^ + qiU)mod s, where 
Pi, Qi e {0, 1, 2, . . . , s — 1}, for all i = 1 to N. Also Cj^ ^ Cj^ and c^j, Cj^ e C implies 
that Cj-^ and Cj^ form the factors of an 0A{N,2, s,2). Since, for all z = 1 to N, Xi 
and Yi are just cyclic permutations of symbols of cj^ and Cj^, the sets {Xj, Fj} 
form the factors of Orthogonal Arrays with parameters (A^, 2,s,2). Therefore, 
Z = [c'l : C2] being a vertical juxtaposition of the matrices [Xi : Yi] for i = 1 to N, is 
an OA{N'^, 2, s, 2). Hence, c[ and 4 form the factors of an OA{N^, 2, s, 2). 

Case 2 : Cj^ — cj^ such that cj^, cj^ 7^ Ck+i 

Xi — {cj + PiU)mod s and Yi — {cj + qiU)mod s, where Cj^ — Cj^ — Cj (say) 
and Pi, Qi G {0, 1, 2, . . . , s — 1}, for alH = 1 to N. Since we know that Cj^ 7^ Ci^, 
and also the fact that and together form the factors of an 0A{N,2, s,2), we 
therefore can conclude that all possible 2-tuples on {0, 1, 2, . . . , s — 1} belong to the 
set {{pi,qi) : i = 1 to N} and also that each tuple occurs exactly N/s times. Hence, 
Z = [c[ : C2] which is a vertical juxtaposition of the matrices [Xi : yj] for i = 1 to N, 
is an OA{N'^, 2, s, 2)and thus c'^ and 4 form the factors of an OA{N^, 2, s, 2). 

Case 3 : Cj^ = Ck+i and Cj^ 7^ Ck+i 

c\ is a juxtaposition of all, a e {0, 1, 2, .... s — 1} with each all appearing N/s 
times, c'2 is a juxtaposition of Cj^ + I3U, /3 e {0, 1, 2, . . . , s — 1} where each Cj^ + 
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has symbols from {0, 1, 2, . . . , s — 1} an equal number of times (since A is orthogonal 
of order 1). Hence c[ and C2 contain all possible 2-tuples over {0, 1, 2, . . . , s — 1} an 
equal number of times. 

Case 4 : Cj^ = Cj^ = Cfe+i 

= implies 7^ Ci^. Also since and contain all possible 2-tuples over 
{0, 1, 2, . . . , s — 1} (since A is orthogonal of order 2), c[ and c'2 contain all possible 
2-tuples over {0,l,2,...,s — 1} equal number of times. 

Since c'j^ and were chosen arbitrarily, C forms the set of all the factors of 
an OA{N^,k'^ + 2k, s, 2). Hence, the construction given in Algorithm 1 gives an 
OA{N'^, k'^ + 2k, s, 2) for every existing OA{N, k, s, 2). 

It can be noticed that the linearity constraint is not necessary for proving that 
the method gives an OA{N'^, k'^ + 2k, s, 2). Hence the above lemma is proved. 

3.4 Generating Non-Linear Orthogonal Arrays 

Although we originally proposed the method for generating linear orthogonal arrays, 
it is possible to generate non-linear orthogonal arrays using the same construction. 
For accomplishing this we must take a linear orthogonal array to begin with, choose 
any one of its columns and cyclically permute its symbols. Then clearly the resulting 
orthogonal array will loose its linearity property. However, even then the construc- 
tion will work well for this kind of seed orthogonal array and generate a non-linear 
orthogonal array. We discuss further about this approach and its correctness. 
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Given a linear orthogonal array A{N, k, s, t), we randomly choose a column q of 
which we shall cyclically permute the symbols. After undergoing a cyclic permutation, 
the new column can be represented as follows: 

c[ — Ci + all 

where o; is a non-zero element of GF{s) and t/ is a unit column vector of dimensions 
X 1. 

Then the columns of the orthogonal array B generated using c[ from the seed 
array A may be of the form: 

c' = {{ci + all) <S)U + U Cj)mods = {ci®U + U ® Cj + all ® U)mods 

Irrespective of whether Cj is a zero column vector or is equal or not equal to c^, 
we always get a symbol-wise cyclic permutation of the column of B we would have 
generated if A were linear. Similarly, the result holds for 

d = ((cj all) + U (S> Cj)mods 

Since the the columns of the array we have obtained is proved to be equal to or a 
cyclic permutation of symbols of the columns of the linear orthogonal array which we 
would otherwise have generated, the array thus produced is orthogonal. Further, since 
it does not contain a zero vector for a row, it is not linear. Hence, we can successfully 
generate non-linear orthogonal arrays using the aforementioned construction. 
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Chapter 4 



Results 



The significance of a new construction lies in its ability to generate new orthogonal 
arrays which can be contributed to the existing libraries. Hence, it is important to 
show that the construction proposed in this work is actually capable of generating new 
orthogonal arrays. Thereby, in this chapter we shall discuss about the new orthogonal 
arrays that can be constructed using the proposed method. Comprehensive libraries 
of orthogonal arrays can be found in Hedayat et al. (1999) or on the following websites: 

http:/ / www2. research, att. com/ njas/ oadir/ 
http://support.sas. com/techsup/technote/ts723.html 

4.1 List of New Contributions 

In this section we shall list out tables of orthogonal arrays which are generated using 
the proposed construction but are not found in any of the above mentioned libraries. 
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Table 4.1: Orthogonal Arrays with Two Levels and Strength 2 



SI. No. 


Seed OA 


Generated OA 


I 


OA(8,5,2,2) 


OA(64,35,2,2) 


2 


OA(8,7,2,2) 


OA(64,63,2,2) 


3 


OA(12,ll,2,2) 


OA(144,143,2,2) 


4 


OA(16,15,2,2) 


OA(256,255,2,2) 


5 


OA(20,19,2,2) 


OA(400,399,2,2) 



Table 4 .2: Orthogonal Arrays with Two Levels and St rength 3 
SI. No. Seed OA Generated OA 

1 OA(24,12,2,3) OA(576,168,2,3) 

2 OA(32, 16,2,3) OA(1024,288,2,3) 

3 OA(40,20,2,3) OA(1600.440.2,3) 

4 OA(48,24,2,3) OA(2304,624,2,3) 

5 OA(56,28,2,3) OA(3136,840,2,3) 

6 OA(64,32,2,3) OA(4096,1088,2,3) 

7 OA(72.36.2.3) OA(5184. 1368.2,3) 



Table 4.3: Orthogonal Arrays with Two Levels and Strength > 3 
SI. No. Seed OA Generated OA 

1 OA(80,6,2,4) OA(6400,48,2,4) 

2 OA(128,9,2,5) OA(16384,99,2,5) 

3 OA(64,7,2,6) OA(4096,63,2,6) 



Table 4.4: Orthogonal Arrays with Three Levels and Strength > 2 



SI. No. 


Seed OA 


Generated OA 


1 


OA(27,13,3,2) 


OA(729,195,3,2) 


2 


OA(81,40,3,2) 


OA(6561, 1680,3,2) 


3 


OA(54,5,3,3) 


OA(2196,35,3,3) 
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Table 4.5: Orthogonal Arrays with More Than Three Levels 



SI. No. 


Seed OA 


Generated OA 


1 


OA(16,5,4,2) 


OA(256,35,4,2) 


2 


OA(64,21,4,2) 


OA(4096,483,4,2) 


3 


OA(64,6,4,3) 


OA(4096,48,4,3) 


4 


OA(25,6,5,2) 


OA(625,48,5,2) 


5 


OA(49,8,7,2) 


OA(2401,80,7,2) 


6 


OA(64,9,8,2) 


OA(4096,99,8,2) 


7 


0A(81, 10,9,2) 


OA(6561, 120,9,2) 
OA(11^168,ll,2) 


8 


0A(121, 12,11,2) 


9 


OA(169, 14,13,2) 


OA(13^224,13,2) 


10 


OA(256, 17,16,2) 


OA(2i6,288,16,2) 
OA(17^360,17,2) 


11 


OA(289,18,17,2) 
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Chapter 5 



Conclusion and Future Work 

In this work, we presented a novel construction algorithm for generating orthogonal 
arrays from seed linear orthogonal arrays. The proposed method works well for seed 
linear orthogonal arrays of all strengths and levels. We also discussed the proof 
of correctness of the construction and also the possibility of generating non-linear 
orthogonal arrays using the same construction. The results show that the proposed 
construction is indeed capable of contributing to the existing libraries of orthogonal 
arrays. Lists of new orthogonal arrays generated using this method have been 
provided in the Results section. 

Extensive experimental observations suggest that the proposed construction works 
well even if the required linearity condition is dropped. This gives us reason to 
conjecture that the construction would work well for any seed orthogonal array. We 
have already proved this result for seed orthogonal arrays of strength 2. Our future 
work will be primarily directed towards proving the correctness of the construction 
for non-linear seed orthogonal arrays of any strength. 
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